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Inelastic x-ray scattering (IXS) with very high (meV) energy resolution has be- 
come a valuable spectroscopic tool, complementing the well established coherent 
inelastic neutron scattering (INS) technique for phonon dispersion investigations. 
In the study of crystalline systems IXS is a viable alternative to INS, especially 
in cases where only small samples are available. Using IXS, we have measured 
the phonon dispersion of Ndi.86Ceo.i4Cu04+<s along the [£,0,0] and [£,£,0] in- 
plane directions. Compared to the undoped parent compound, the two highest 
longitudinal optical (LO) phonon branches are shifted to lower energies because of 
Coulomb-screening effects brought about by the doped charge carriers. An addi- 
tional anomalous softening of the highest branch is observed around q = (0.2,0,0). 
This anomalous softening, akin to what has been observed in other compounds, 
provides evidence for a strong electron-phonon coupling in the electron-doped high- 
temperature superconductors. 

1 Introduction 



Inelastic neutron scattering (INS) measurements of phonon dispersions inJjhc p- 
type superconductors Lai^Sro.isCuO^atl&tM, oxygen-doped La2CuC>4 +< 5la, and 
YBa2Cu3O6+50&0 show an anomalous softening of the highest longitudinal opti- 
cal (LO) phonon branch with doping, especially along q = [£,0,0] in these p-type 
cuprate superconductors. The atomic displacement assigned to this branch is the 
Cu-0 bond-stretching (BS) modecra. The. interpretation of such an anomaly in- 
vokes a strong electron-phonon couplingcra, but—itJa-still a challenge to provide a 
quantitative model of this coupling (see Ref. BBmJ'EEl). 

Similar anomalies also develop with doping in other perovskite oxides, 
such as Lai_ x Sr x Mn03Lil and Bai_ x K x Bi03E3, which are nearly cubic, and in 
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La2_ x Sr x NiOji3, which is iso-structural to Lai.ssSro.isCuC^+s. In the case of 
Bai_ x K x Bi03, there is evidence that the strong electron-phonon coupling is re- 
sponsible for both the anharmonic behaviour of the high-energy LO modes and 
the high superconducting transition temperature (30 K) of this materialtij. How- 
ever, a microscopic model capable of explaining the electronic and lattice dynamical 
properties of Bai_ x K x BiC>3 is still lackingcj. 

The effect of an anomalously strong electron-phonon coupling is still unclear. 
At this point, one may not even exclude the possibility that it may be pair-breaking 
for the d-wave superconducting order parameteiO. Therefore, it is very important 
to analyse the phonon anomalies in as many cuprate families as possible, and to 
compare these anomalies with the electronic properties. In this context, the dif- 
ferent character of electron and-hole carriers is of central importance: Cu 3d x 2_ y 2 
(O 2p) foiij|i(p)-type cupratestl This difference leads to very different electronic 
structureaiH. Since the phonon anomalies are related to a coupling between charge 
fluctuations and phonons, charge carriers with different character may induce quite 
different electron-phonon interactions. Therefore, it is imperative to extend studies 
of the electron-phonon coupling to the n-type cuprates. _ . 

X-ray and neutron scattering are complementary experimental toolalj. An ad- 
vantage of IXS over INS is the ability to work with very small samples due to the 
very small beam size and high brilliance obtainable at 3 rd generation synchrotron 
sources. Moreover, the scattering volume for x-rays in the energy range of 10 to 20 
keV is strongly limited by photoelectric absorption for Z > 4. In the case of the 
cuprates, IXS probes a depth between 10 and 100 /im, comparable with the lateral 
size of the x-ray beam obtained on the ID16 (or ID28) beam line at the ESRF. This 
allows experiments with scattering volumes as small as 10~ 4 mm 3 . Therefore, IXS 
is the only technique available for phonon dispersion measurements of small sam- 
ples. This is relevant for the electron-doped superconductor Nd2- x Ce x Cu04 + 5 in 
which the oxygen non-stoichiometry S may vary throughout the sample as a con- 
sequence of the reduction procedure required to create the superconducting phase. 
The small scattering volume ensures that effects due to any composition gradient 
are minimised. A further advantage of the short photoelectric absorption length is 
the reduction of multiple scattering processes^ 

Another important difference between INS and IXS is that the magnetic cross 
section is negligible for x-rays, whereas it is comparable to the nuclear cross section 
for neutrons. Therefore, x-rays are essentially insensitive to magnetic excitations in 
the THz regime, which can be helpful for separating magnetic and lattice contribu- 
tions. Indeed, in Nd2- x Ce x Cu04+a there are very strong crystal field excitations 
with a neutron cross section about ten— times greater than that for phonons, in 
particular in the range from 4 to 8 THzll3. 

Other important differences between IXS and INS are: a) The x-ray cross sec- 
tion is highly coherent. In particular, the entire elastic contribution comes from 
<S(AQ, Az/), where AQ and Av are given by the instrumental resolution. There- 
fore only static or quasi-static structural and/or chemical disorder contributes to 
the quasi-elastic signal. For perfect, pure crystals no elastic signal is present far 
from the Bragg reflections; b) AQ and Av are decoupled over a wide range of 
(Q, f). Moreover, the frequency (energy) Av is essentially independent of the in- 
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cident frequency (energy) z/,-. This allows access to high Q and/or high frequency 
transfers (> 20 THz) with good resolution in both v and Q. Because of the highly 
collimatcd incident beam, high Q-resolution is easily obtained; c) the shape of the 
IXS instrumental energy resolution is not Gaussian as it is for a neutron triple-axis 
spectrometer, but Lorentzian. This strongly enhances the tails of both the elastic 
and inelastic signal, and tends to prohibit IXS from measuring phonons in the pres- 
ence of a large elastic signal such as Bragg peaks or strong diffuse scattering. As a 
consequence, inelastic data are usually not available at the zone centre; d) in the 
case of high-energy optical branches in the cuprates, the relevant scattering is due 
to oxygen motion, while the relevant absorption cross section is due to neodymium. 
As a consequence, the signal from the high-energy optical branches, which is rela- 
tively low because of the ~ terms of the dynamical structure factor and of the Bose 
term, is further depressed compared to the signal of the low-energy modes and the 
elastic signal, to which all the atoms contribute. The strongly enhanced tails of the 
elastic and low-energy inelastic signal (point c)) and Z contrast (point d)) above 
have posed great difficulties for IXS measurements of high-energy optical branches 
in transition metal oxides. We will show in the next section that we have been 
able to overcome these difficulties through a careful choice of thermodynamic and 
kinematic conditions. 

2 Experiment and Results 

The inelastic x-ray scattering experiment was carried out at the very-high- 
energy-resolution beam-line ID16 at the European Synchrotron Radiation Facility 
(ESRF). X-rays from an undulator source were monochromated using a Si (111) 
double-crystal mpnochromator, followed by a high-energy-resolution backscattering 
monochromator Ej operating at 15816 eV (Si (888) reflection order) and 17794 eV 
(Si (999) reflection order). A toroidal gold-coated mirror refocused the x-ray beam 
onto the sample, where a beam size of 250 x 250 fan 2 full-width-half-maximum 
(FWHM) was obtained. The scattered photons were energy-analysed by a spheri- 
cal silicon crystal analyser of 3 m radius, operating at the same Bragg reflection as 
the monochromator EJ. The total energy resolution was 1.6 THz (6.6 meV) FWHM 
for the Si (888) reflection and 1.1 THz (4.4 meV) FWHM for the Si (999) reflec- 
tion. The momentum transfer Q was selected by rotating the 3 m spectrometer 
arm in the scattering plane perpendicular to the linear x-ray polarisation vector 
of the incident beam. For most of the scans, the momentum resolution was set 
to ps 0.087 A -1 in both the horizontal and the vertical direction by an aperture 
of h X v = 20 X 20 mm 2 of the slits in front of the analyser. In order to achieve 
higher Q-resolution, for some selected scans we set the aperture to 20 x 10 nim 2 . 
and 10 x 10 mm 2 . Further experimental details are given elsewhere (see Ref. E2lEil 
and references therein). 

The samples are two single crystals grown by the travelling-solvent floating-zone 
method in 4 atm of O2 at Stanford University: a crystal of the parent insulator 
Nd2Cu04+,5 and a crystal of the Ce-doped superconductor Ndi.geCeo.uGuO^^. 
The latter sample was oxygen-reduced after the growth and had a superconducting 
onset temperature of T c = 24 K. Further details are given elsewherec3. The samples 
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are of very good crystalline quality. Initial characterisation work, performed on 
beam- line 7-2 at the Stanford Synchrotron Radiation Laboratory (SSRL), gave 
rocking curve widths of 0.02° (FWHM) around the [h, 0,0] direction. 

The IXS experiment was performed in reflection geometry, and the probed scat- 
tering volume corresponded to about 1.5 x 10~ 3 mm 3 . IXS scans were performed 
in the -2< v <24 THz range, in the r = (6,0,0), r = (7,1,0), and r = (5,5,0) 
Brillouin zones. ^| The samples were mounted on the cold finger of a closed-cycle 
helium cryostat, and cooled to 15 K. The low temperature and high momentum 
transfer were chosen so as to optimise the count rate on the high-frequency optical 
mode while limiting the loss of contrast due to the contribution from the tails of 
the intense low-frequency acoustic modes. The effect of lowering the temperature 
to 15 K and of the high momentum transfer is to enhance all the signals by the 
Debye-Waller factor, and to reduce the intensities of the low-energy (kT < hv < 
kTfl-r) phonons by the Bose factor. This is shown in Figs. [I] (a) and (b), which 
compare data obtained at 300 K and 15 K. We note that the main additional con- 
tribution at high energy transfers comes from the tail of the longitudinal acoustic 
(LA) mode (dashed line) , which is strongly reduced at low temperature (see Fig. [I] 
(b))- 

The two highest-energy modes for Ndi.g6Ceo.i4Cu04 + 5 and Nd2Cu04+5 at 
15 K and at the reduced wave-vector q = (0.1,0.1,0) are shown in Fig. |l| (b) 
and (c), respectively. We observe a shift towards lower energy from the parent 
antiferromagnetic insulator to the Ce-doped superconductor. The two highest LD 
branches are assigned to the Cu-0 bond-stretching and 0(2) vibration modesotl, 
respectively. While this assignment is correct near the zone centre, its validity over 
the whole Brillouin zone is discussed in Sec. ||. However, for the sake of simplicity, 
we refer to the branches in this way in the text below. 

In Fig. |], we also show the low-energy acoustic and optical longitudinal modes 
at the zone-boundary (M point) in N^CuC^+a (d) and near the zone-boundary 
(q = (0.45,0.45,0)) in Ndi.86Ceo.i4Cu04+«5 (e) and (f). In the doped specimen, 
measurements at the M-point were not possible because of the strong elastic signal 
at t = (5.5, 5.5, 0) attributed to a secondary structural phase formed during the 
reduction procedure. Near the zone boundary, the strongest mode is expected to be 
the lowest-energy longitudinal optical mode. We observe the lowest optical mode 
near the zone-boundary around 6 THz in both undoped Nd2Cu04+,5 and doped 
Ndi.86Ceo.i4Cu04+5. This observation is different from that previously reported 
for Nd2Cu04+5 based on INSEl, where the same mode was found at lower frequency 
(about 4 THz). We note that-in this frequency range, a very strong crystal field 
excitation is observed by INSE3. In order to confirm our assignment we performed 
a scan with higher energy resolution, which is shown in Fig. |l| (f)^|. 

The frequencies of the peak positions extracted from these and many other 
scans along the main in-plane symmetry directions (A = (£,0,0) and £ = (£,£,0)) 
are summarised in Fig. |^. The IXS results for Ndi.86Ceo.i4Cu04+5 (°) and 
Nd2Cu04+5 (•) are shown together with the INS data for the highest Cu-0 BS 



a All reciprocal lattice vectors are expressed in reciprocal lattice units. 

b Note that we observe three acoustic contributions (both longitudinal and transverse) at low 
frequency due to the finite Q-rcsolution of our experiment. 
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Figure 1. IXS phonon spectra. Left panel: for q = (0.1,0.1,0) in Nd1.86Ceo.14CuO4.1_ g (o) at (a) 
300 K and (b) 15 K and (c) Nd 2 Cu0 4+ a (•) at 15 K, with a resolution of Av = 1.6 THz. Right 
panel: (d) at the reduced wave-vector q = (0.5, 0.5, 0) in Nd2Cu04+a (•) at 15 K, with a resolution 
of Av =1.6 THz, and at the reduced wave-vector q = (0.45,0.45,0) in Ndi.ggCeo.i4Cu04+i (o) 
at 15 K with a resolution of (e) Av =1.6 THz and (f) Av = 1.1 THz. Lines indicate the best fit 
with a harmonic oscillator model. The dashed lines in the left panels indicate the contribution 
from the tail of the LA mode. 



mode (A)H. The lines are the results of a numerical calculationE-l based on a shell 
model, where the inter-atomic potentials are chosen following Ref. c3. We simu- 
lated the phonon modes in the parent antiferromagnetic insulator Nd 2 Cu04 + a by 
adding an unscreened Coulomb interaction (dot-dashed lines). The results repro- 
duce very well thje experimental dispersion of the undoped parent compound as 
measured hy.INSI-J and IXS. For the comparison with Ndi.86Ceo.i4Cu04+,5, follow- 
ing Ref. E3, we added a screened Coulomb interaction (solid lines) in order to 
simulate the effect of the free carriers introduced by dopingc3 . 



3 Discussion 



The zone centre. — The shift at the zone centre of the high-energy phonon 
branches of Ndi.seCeo.uCuC^+awith. respect to Nd2Cu04 + <sis due to the closing 
of a large longitudinal optical - transverse optical (LO-TO) splitting (also known 
as Lyddane-Sachs-Teller splitting), which is very well reproduced by our numer- 



sattll 'dastuto'cm: submitted to World Scientific on February 1, 2008 5 



£ A 




0.5 1 

B,o,o] 

Figure 2. Dispersion of the longitudinal phonon modes in Ndi.86Ceo.i4Cu04+a (o) and 
Nd2Cu04+a (•) measured by IXS. For comparison, we include the dispersion of the highest LO 
mode in Nd2Cu04_|_<5 (A) as measured by INS (from Ref. 1). Solid (dot-dashed) lines indicate 
the results of a lattice dynamics calculation with a screened (unscreened) Coulomb interaction. 
The shaded areas indicate regions where the dynamic structure factor is very weak and the details 
of the second highest mode can not be determined. 



ical calculations. The corresponding Ai and A3 (or Si and S3) branches in 
Nd2Cu04+5 are separated by almost 3 THz at the zone centre, as observed by INStl, 
while for undoped La2Cu04+5c3 the LO and TO branches are very close. Hence, 
there does not exist a similarly strong screening effect in Lai.ggSro.isCuO^a for the 
high-energy in-plane LO modeso. The screening vector k s , which we found to be ap- 
proximately 0.39 A , is comparable to the one estimated for Lai.ssSro.isCuO^aEa. 

Although our calculation with a screened Coulomb potential reproduces the 
closing of the LO-TO splitting, we still observe an anomalous additional softening of 
the highest LO branch near q = (0.2, 0, 0) and above approximately q = (0.2, 0.2, 0) 
which is not reproduced by our calculations (see Fig. |^). 

A (or [£,0,0],) line. — Along the A line ([£,0,0]; right side of Fig. g), 
the highest LO branch softens in frequency by about v\ — V2 ~ 1.5 THz between 
qi = (0.1,0,0) and q2 = (0.2,0,0), which is a shift cotmiiarable to the anomalous 
shift observed in Lai.ssSrn.isCuO^a at q = 0.25-0.3 aOa Therefore, we believe 
that this anomalous softening is of the same nature as the one observed in p-type 
Lai.85Sro.i5Cu04+,5. For q > 0.2, the highest A branch recovers rapidly. This be- 
haviour is most likely due to an anti-crossing with the second highest branch which 
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is mainly associated with vibrations of the 0(2) atoms in the ^-direction. There- 
fore, within a standard anti-crossing framework, one should interpret the highest 
longitudinal intensities for q = (0.3,0,0) and above as being mainly due to 0(2) 
vibration. Within this scenario, the second-highest branch increases in frequency 
in the middle of the zone as in the un-doped compound. Thus, except for the fact 
that the LO-TO gap closes, this mode seems to be insensitive to doping. The LO 
bond-stretching mode just above q = (0.2,0,0) is then found at rather low ener- 
gies, ~12 THz, with a flat-dispersion, as has been reported for Lai.85Sro.isCu04+5 
for q > 0.25a or q > 0.3a, but can not be unambiguously followed to larger q- 
values. However, the measured intensities do not behave as expected for a simple 
anti-crossing picture (see Ref. E3). In an anti-crossing scenario, one expects an 
exchange of character of the two branches, which would imply an exchange in the 
intensities. Instead, we observe an intensity that agrees with our calculation, which 
does not include any coupling between the modes, and hence does not allow for an 
anti-crossing. This appears to reflect the fact that we still lack a realistic model 
of the phonon anomaly. Furthermore, in the (Q,w) region indicated by the shaded 
area, the dynamic structure factor is very weak and the mode of interest is closer to 
the tails of the low-energy.-|modes. Therefore, we canjjjOt, resolve broadening effects 
as in the hole-doped casecra, or even band splittingc3 : tL3. Consequently, with the 
present quality of IXS data, the position for the second-highest branch can only be 
defined as the centre-of-mass frequency in this (Q,w) region. 

S (or [f, £,0])) line. — Along the £ line (or [f, £, 0]; left side of Fig. |), the 
softening of the highest LO branch between q = (0.1,0.1,0) and q = (0.2,0.2,0) 
is much less than that observed at q = (0.2,0,0). However, here we observe a 
softening of the second branch, down to 13 THz for £ > 0.2. In the framework 
of the anti-crossing picture, this has to be interpreted as the softening of the Cu-0 
BS mode, but the same caution must be taken in the analysis of this lower-energy 
mode in the (Q,w) region indicated by the shaded area as in the analysis of the 
mode along A. However, we note that: i) due to the different symmetry of the 
orbitals involved in n- and p-type doping, it has been predicted that ±he bond- 
stretching mode would be particularly soft at the [£, £, 0] zone boundary^, and ii) 
the second phase observed in the superconducting sample (see Sec. ||), with peaks 
at r = (h + 5, h + i, 0), may influence the dynamics near the zone boundary. 

4 Conclusions 

The present results on the in-plane longitudinal phonon dispersion in 
Ndi.86Ceo.i4Cu04 + 5, when compared to undoped Nd2Cu04 +( 5EJ, show that the 
anomalous softeninft.af-.the highest LO phonon mode previously observed in hole- 
doped compoundsclufflla-El is also pxesfijxLin the electron-doped cuprates. This 
observation supports the hypothesisaffla-BO of an electron-phonon coupling origin 
of this feature. 

These results furthermore demonstrate that high-energy resolution inelastic x- 
ray scattering has developed into an invaluable tool for the study of the lattice 
dynamics of complex transition metal oxides, allowing measurements on small single 
crystals which can not be studied by the traditional method of inelastic neutron 
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scattering. 
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